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Abstract. We show that the inclusion of special relativistic corrections in the revised OPAL and MHD equations 
of state has a significant impact on the helioseismic determination of the solar age. Models with relativistic 
corrections included lead to a reduction of about 0.05 — 0.08 Gyr with respect to those obtained with the old 
OPAL or MHD EOS. Our best-fit value is t so i s = (4.57 ± 0.11) Gyr which is in remarkably good agreement with 
the meteoritic value for the solar age. We argue that the inclusion of relativistic corrections is important for 
probing the evolutionary state of a star by means of the small frequency separations Sue, n = f<,n — t^+a.n-ij for 
spherical harmonic degrees £ — 0, 1 and radial order n>!. 
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1. Introduction 

The possibility of using helioseismology to constrain the 
solar age has been discussed by several authors in the 
past. Very recently Dziembowski et al. (1999) have shown 
that the most robust and accurate method is provided by 
the small frequency separation analysis (SFSA), (5^y.„ = 
vi,n~ ^+2.71-1 , for spherical harmonic degrees i = 0, 1 and 
radial order n > £ (Tassoul 1980). 

The important property of this quantity is its strong 
sensitivity to the sound-speed gradient near the solar cen- 
tre and its weak dependence on the details of the treat- 
ment of the outer layers. Despite our ignorance of a reliable 
convection model for the solar envelope we are therefore 
able to verify how well our models are able to reproduce 
the deep radiative regions, in particular the solar core. 
Since the properties of the core are mainly determined by 
the present central hydrogen abundance, and the latter 
is influenced by the solar age, SFSA is a reliable tool to 
examine the seismic age of the Sun. 

Adopting the OPAL equation of state ( Rogers et al. 
1996) a seismic age of (4.66±0.11) Gyr has been obtained 
by Dziembowski et al. (1999), which is consistent with the 
meteoritic age (4.57 ± 0.02) Gyr of Bahcall et al. (1995). 

The aim of this paper is to show that an important in- 
gredient in this type of analysis is the usage of an accurate 
equation of state (EOS). In particular, by the inclusion of 
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the special relativistic corrections, like in the updated ver- 
sion of the OPAL EOS, the helioseismic age of the Sun is 
reduced to (4.57±0.11) Gyr, which is in remarkable agree- 
ment with the meteoritic value. 

Elliott & Kosovichev (1998) have demonstrated that 
the inclusion of relativistic corrections in the EOS leads 
to a better agreement between the solar models and the 
seismic Sun. By inverting SOI-MDI/SOHO p-mode fre- 
quencies they found that the solar adiabatic exponent 
Ti is much better reproduced by solar models including 
the relativistic contribution to the Fermi-Dirac statistics. 
Since the improved EOS causes a decrease of 0.2% in the 
adiabatic index T% in the solar centre, the sound speed 
(oc yTT) is reduced by about 0.1%. Therefore, the influ- 
ence of the relativistic corrections should also be visible 
in the small frequency separations bv^ n . Indeed, Bonanno 
et al. (2001) have found that including this effect in the 
value of Ti improves the agreement in bv^ n between solar 
models and observations, thereby confirming the results 
of Elliott & Kosovichev (1998). 

In addition to the age, the central hydrogen abun- 
dance is also crucially dependent on the precise value 
of S pp (0), the zero-energy astrophysical S-factor for the 
proton-proton fusion cross section. Schlattl et al. (1999) 
and Antia & Chitre (1999) have shown, using the old ver- 
sion of the OPAL EOS, that an increase of S pp (0) by about 
4% with respect to Adelberger et al.'s (1998) value yields 
a better agreement with the observed frequencies for an 



2 



A. Bonanno et aL: The solar age and the relativistic correction 



Table 1. Characteristic quantities of selected solar models. The indices 0, ph, cz, and c denote initial, photospheric, 
bottom of convective envelope, and centre, respectively. MHD-R is the abbreviation for the MHD EOS containing the 
relativistic corrections in I\. 



Model 


age 


EOS 


Y 


Z 






r cz 


x c 


Y c 


Pa 


T c 


ti pp {0) 


Gyr 




gem - 3 


10° K 


10- 25 McVb 


1 


4.58 


OPAL 01 


0.2755 


0.01995 


0.2453 


0.01805 


0.7132 


0.3353 


0.6432 


152.87 


15.73 


3.89 


2 


4.58 


OPAL 96 


0.2749 


0.01995 


0.2449 


0.01806 


0.7132 


0.3289 


0.6428 


152.70 


15.72 


3.89 


3 


4.60 


OPAL 01 


0.2752 


0.01995 


0.2451 


0.01805 


0.7125 


0.3342 


0.6443 


153.16 


15.73 


3.89 


4 


4.60 


MHD-R 


0.2757 


0.01997 


0.2452 


0.01805 


0.7141 


0.3341 


0.6444 


153.22 


15.74 


3.89 


5 


5.00 


OPAL 01 


0.2714 


0.02013 


0.2405 


0.01816 


0.7082 


0.3133 


0.6650 


159.82 


15.84 


3.89 


6 


4.58 


OPAL 01 


0.2758 


0.01989 


0.2460 


0.01803 


0.7118 


0.3362 


0.6423 


151.35 


15.66 


4.00 



age of 4.57 Gyr. For this reason we consider in our analysis 
also different values of S pp (0). 

Including the updated OPAL EOS the best agreement 
between meteoritic and seismic age could be achieved with 
Adelberger et al.'s (1998) S pp {0) = 4.00 x 10~ 25 MeVb. 
Hence, by taking into account the relativistic corrections 
in the EOS there is no need for an artificial increase of 
S pp (0), as suggested by previous works, in order to obtain 
a better agreement between seismic and meteoritic age. 

The code and physics used to compute the various solar 
models are described briefly in the next section, followed 
by the consequences for the seismic age obtained by means 
of the SFSA (Section 3). In the final part the results are 
discussed. 

2. The new solar models 

We computed a large number of solar models using the 
GARching SOlar Model (GARSOM) code which has been 
described in its latest version in Schlattl (2001). Our stan- 
dard model has been compared with other contemporary 
solar models by Turck-Chieze et al. (1998), who found a 
good agreement between various programs. 

The solar photospheric radius and luminosity have 
been assumed to be 695.51 Mm (Brown & Christensen- 
Dalsgaard 1998) and 3.86 46 x 10 33 erg/s, respectively. The 
surface metal ratio has been taken from Grevesse & Noels 

(1993) , thus Z/X = 0.0245. The mixing length parame- 
ter (Bohm-Vitense 1958), initial helium and metal content 
have been adjusted in all models to reproduce these values 
with an accuracy better than 10~ 4 . 

In the actual calculations the latest OPAL-opacities 
(Iglesias & Rogers 1996) completed in the low- 
temperature regime by tables of Alexander & Fergusson 

(1994) have been implemented. The outer boundary con- 
dition was determined assuming an Eddington grey atmo- 
sphere. Microscopic diffusion of hydrogen, helium and all 
major metals is taken into account. For the EOS we used 
either the OPAL- (Rogers et al. 1996) or the MHD-tables 
(Hummer & Mihalas 1988, Mihalas et al. 1988, Dappen 
et al. 1988). The original OPAL EOS (OPAL96) has been 
updated by treating electrons relativistically and by im- 
proving the activity expansion method for repulsive inter- 
actions (Rogers 2001), denoted OPAL01 in the following. 



In the case of MHD EOS the relativistic corrections 
are not directly included in the tables. We have therefore 
corrected the adiabatic index Ti employing the expression 
of Elliott & Kosovichev (1998), 

5Ti _ r^d-rx _ 2 + 2X kT 

Ti ~ ri 3 + 5X m e c 2 ' [ ' 

where T is the temperature, m e the electron mass, c the 
light speed in vacuum, k the Boltzmann constant, and X 
the hydrogen mass fraction. As expected, the correction to 
Fi is negative, since its value is 5/3 for the non-relativistic 
and 4/3 for the extremely relativistic case. 

The nuclear reaction rates are taken either from 
Bahcall et al. (1995) or from Adelberger et al. (1998) 
with S pp (0) being 3.89 x 10~ 25 MeVb in the first and 
4.00 x 10~ 25 MeVb in the latter case. Other differences in 
the reaction rates are not very significant in determining 
the evolutionary stage of the solar core. 

3. Results for the solar age 

We have computed solar models following the evolution 
from the zero-age main sequence with ages ranging from 
4.40 to 5.00 Gyr in steps of 0.1 Gyr. Some basic quantities 
of a selection of models are summarized in Table [j]. 

For the higher ages the initial helium content has to be 
reduced to obtain the correct solar luminosity (compare 
models 1 and 5). Nevertheless, a larger lifetime leads to a 
steeper He profile toward the centre causing a larger cen- 
tral He abundance. The consequent increase of the opacity 
near the core demands an higher central temperature to 
produce the same amount of energy. This effect is fur- 
ther enhanced by diffusion which is operating longer for 
greater ages and is further increasing the central He con- 
tent. Since the relativistic correction to Ti increases with 
temperature (Eq. [I]), the inclusion of relativistic effects 
has a larger influence on older models. The relative differ- 
ences in the profiles of Ti and the density are shown in 
Fig. | 

Models with greater S PP (Q), but the same age, have a 
smaller T c (see models 1 and 6 in Table Q), as the hydrogen 
burning in the core is more efficient. 

In order to determine the seismic age, we calculated 
for all the solar models the small frequency separations 
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Fig. 1. The differences of Ti (solid line) and p (dashed) 
between two models which either neglect (model 1) or 
contain the relativistic corrections (model 2) in the sense 
(model 2 — model 1)/ (model 1). 

Table 2. The best- fit age and the corresponding minimum 
of x 2 f° r the grid with different equations of state and 
different values of S pp (0) in units of 10~ 25 MeVb. 



EOS 


Spp(0) 


I = 

^scis 


xl 


i = 1 

^scis 


xi 


OPAL96 


3.89 


4.664 ± 0.088 


1.05 


4.672 ± 0.088 


1.66 


OPAL01 


3.89 


4.584 ± 0.088 


1.45 


4.624 ± 0.072 


1.66 


MHD 


3.89 


4.664 ± 0.080 


1.00 


4.680 ± 0.095 


1.65 


MHD-R 


3.89 


4.608 ± 0.040 


1.07 


4.640 ± 0.088 


1.25 


OPAL01 


4.00 


4.552 ± 0.080 


1.34 


4.584 ± 0.080 


1.47 



8vi, n for £ = 0,1 and n ^S> £■ These values have been 
compared with latest GOLF/SOHO data for I = 0, 1, 2, 3, 
which have been obtained from long time series, and where 
the asymmetric line profile has been taken into account 
during the data reduction (Thiery et al. 2000). Only the 
frequencies of the mean multiplet (m=0) are used, as for 
them the influence of rotation is smallest. 

For the analysis, the x 2 method has been used, as in 
Dziembowski et al. (1999) or Schlattl et al. (1999); 



1 



M - 



A I 

E 



(<^,n,© — Svi,n, model) 2 



(2) 



with M = 31 for £ = and M = 27 for 1 = 1, and m being 
10 in both cases. It is interesting to notice that including 
the relativistic corrections leads to a reduction of 6i>g n of 
about 0.1 /iHz for low frequencies (Fig.^). 

The results for the % 2 -values in models with different 
ages are shown in Figs. and |]. The best-fit age given by 
the minimal x 2_ value (Xm; n ) and the error determined by 
the condition \ 2 ~ Xmin — 1 are summarized in Table ^[ 

Regardless of whether MHD or OPAL EOS is used, the 
best-fit age is reduced by about 0.05-0.08 Gyr when the 
relativistic corrections are included. The minimal value of 
Xq is not significantly different for all the cases, although 



Fig. 2. The differences of the quantity 8v n ,i between 
two models which either neglect (model 1) or con- 
tain the relativistic correction (model 2) in the sense 
(model 2 — model 1) for an age of 4.20 Gyr (solid line) and 
4.70 Gyr (dashed line). 
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Fig. 3. Xo f° r models with different age, neglecting 
(dashed line) or including the relativistic correction (solid 
line). The MHD-EOS has been used for the models in the 
left panel, the OPAL-EOS in the ones of the right panel. 



the models with OPAL96 EOS have a slightly smaller 
Xo,min than those obtained with OPAL01 EOS. 



It is worth noticing that with 5^(0) = 4.00 x 
10 -25 MeVb the minimum x 2 -value slightly improves 
for both £ = and £ = 1 (Table |). Using OPAL01 
EOS, which includes the relativistic corrections in a con- 
sistent way, we obtain in this case as the best-fit age 
^scis = (4.57 ± 0.11) Gyr, where we have taken the mean 
of the best-fit value for £ = and I — 1 . This provides our 
most reliable value for the seismic solar age. 

4. Conclusions 

By using updated versions of the OPAL and MHD EOS 
the seismic age of the Sun has been redetermined using 
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Fig. 4. Same as Fig. |[ but for Xv 



SFSA with the latest GOLF/SOHO data. The important 
new ingredient in both equations of state is the inclusion 
of the special relativistic corrections. In both cases almost 
the same age has been obtained. 

A crucial quantity in the determination of the seismic 
age is the proton-proton fusion rate. With the older ver- 
sions of the equations of state, a rate about 4% higher 
as the value of Adelberger et al. (1998) appears to be 
favoured, in order to obtain a better agreement between 
seismic and meteoritic ages. However, with the updated 
versions of the OPAL and MHD EOS the seismic age ob- 
tained with Adelberger et al.'s (1998) value for S pp (0) is 
(4.57±0.11) Gyr, which is in excellent agreement with the 
meteoritic age of 4.57 Gyr (Bahcall et al. 1995). 

Therefore, the presently favoured value for 5*^^(0) is 
4.00 x 10~ 25 MeVb. However, since the uncertainties, in 
particular, in the opacities are supposed to be of the or- 
der of a few percent, S pp (0) can only be determined with 
a similar accuracy by comparing seismic and meteoritic 
ages. 

A further source of uncertainty is the centrifugal and 
magnetic distortion, but these effects can be neglected for 
the Sun, as discussed by Dziembowski et al. (1999). 

We expect to have asteroseismic data on solar-type 
stars with a precision of about 0.1 /iHz from future space 
missions or high-precision ground-based multi-site spec- 
trographs observations. We thus think that this effect 
must be included in the standard modelling of solar-like 
stars when discussing the evolutionary changes in the stel- 
lar core. 
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